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Abstract 

In this paper, the existence and uniqueness of weak solution for a generalized Boussinesq 
equation that couples the mass and heat flows in a viscous incompressible fluid is considered. The 
viscosity and the heat conductivity are assumed to depend on the temperature. The boundary 
condition on velocity of the fluid is non-standard where the dynamical pressure is given on some 
part of the boundary, and the temperature of the fluid is represented in a mixed boundary 
condition. 
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1 Introduction 

The existence and uniqueness of the solution for the generalized Boussinesq equation have been 
studied extensively by many authors, see O [6l [HI [9l [TO] and the references therein. However, the 
boundary conditions in these works are Dirichlet type for velocity of the fluid. In [4] the blow-up 
and global existence for nonlinear parabolic equations with Neumann boundary conditions is con- 
sidered. The recent works on the existence and uniqueness of the solution of Boussinesq equation 
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Program of China (2011CB808002), and the National Research Foundation of South Africa. 
^The corresponding hor. Email: bzguo@iss.ac.cn 



can be found in [3l O [151 US]- [S] studies the existence and uniqueness of stationary Boussinesq 
system with non-smooth mixed boundary conditions for the temperature, and non-smooth Dirich- 
let boundary condition for the velocity. The local existence and uniqueness of the solution for the 
boundary-value problem for the stationary Boussinesq heat and mass transfer equations with the 
inhomogeneous non-standard boundary conditions for the velocity, and the mixed boundary condi- 
tions for the temperature and concentration are obtained in |15j . [7] investigates the existence and 
uniqueness of weak solutions of the stationary Boussinesq system with the homogeneous Dirich- 
let boundary conditions for the velocity and temperature. The the existence and uniqueness of 
weak solutions of the time-periodic solutions for a generalized Boussinesq equation with Neumann 
boundary conditions for temperature are studied. The global regularity of the classical solutions 
for a 2D Boussinesq equation with the vertical viscosity and diffusivity has been investigated. 

In this paper, we are concerned with the existence and uniqueness of the solution of a gen- 
eralized Boussinesq equation with nonlinear diffusion of the velocity and temperature, where the 
press boundary condition on the velocity of the fluid and the mixed boundary condition on the 
temperature are given, The system is described by the following initial-boundary conditions: 

dz{x, t) 



dt 



j{w{x,t))Az{x,t) + {z{x,t),V)z{x,t) 

= fl{x,t)-VTT{x,t) 

div(z(x,t)) = 
dw{x, t) 



I3gw{x,t) 



dt 



div{k{w{x,t))Vw{x,t)) + {z,V)w = f2{x,t) 



z{x,t)r = 0, 7r(x,t) + = w{x,t) = 



z{x, t) = 0, 



uvoix t] 

-k{w{x,t)) ^ ' ^ =V2{x,t) 

O'" 



n 



z{x,0) = zi{x), ?i;(x,0) = Wq{x) 



on Q, 
on Q, 

on Q, 

on Si, 

on S2, 
on 0, 



(1.1) 



where C {N = 2,3) is a bounded domain with the smooth boundary F = Fi U r2 where 
Fi n F2 7^ 0, Q = 1^ X (0,T) for the given T > 0, Sj = Fj x {0,T),i = 1,2, S = F x (0,r), n is 
the normal vector exterior to F, z{x, t) G denotes the velocity of the fluid at x € O and time 
t € (0,T), TT{x,t) € M is the hydrostatic pressure, w{x,t) € M is the temperature, fi{x,t) is the 
external force, g is the gravitational vector function, 7(-) > is the kinematic viscosity, k(-) > is 
the thermal conductivity, /3 is a positive constant associated to the coefficient of volume expansion, 
f2{x,t) is the heat source strength, and vi,V2 are prescribed functions to be determined later. The 
i-th component of {z{x,t),Vz{x,t))z(x,t) in Cartesian coordinate is given by 



N 



{{z{x,t),V)z{x,t))i 



Y^Zj{x,t)^^4z^^ {z{x,t),V)w{x,t) 



dxj 



N 

^Zj{x,t] 



dw{x, t) 
dxi 



and 



Zt{x, t) = z{x, t) — Zn{x, t) ■ n, Zn{x, t) = {z{x, t) ■ Tl). 
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The classical Boussinesq equation where 7 and k are constants has been studied widely in literature, 
see for instance [HI [12] . Equation p.ip that is paid less attention represents the physical system for 
which we the variation of the fluid viscosity, and the thermal conductivity with temperature cannot 
be ignored (see e.g., [10] and the references therein). Some researches on the generalized Boussinesq 
equation and the Boussinesq system with nonlinear thermal diffusion are available in [SJ [6l [HI [9] but 
the boundary conditions in this works are homogeneous. It should be pointed out that the variation 
of the viscosity with the temperature is important in understanding the process of the flow. As 
a first step to this point, we establish, in this paper, the the existence of non stationary solution 
of Equation (|l.ip where the dynamical pressure is given on some part of the boundary, and the 
boundary condition for temperature of the fluid is mixed. Due to the stronger nonlinear coupling 
between the equations. Equation p.ip is much difficult than the classical Boussinesq equations. 

We proceed as follows. In section 2, we give some preliminary results on Equation p.ip . With 
the Galerkin method, we show, in section 3 the existence of the global solution to the problem 
(jl.ip under certain conditions on the temperature dependency of the viscosity and the thermal 
conductivity. The assumption on non-constant gravitational field would be helpful in some other 
geophysical models. The uniqueness of the weak solution is presented in section 4. 



2 Preliminary results 



For notation simplicity, we do not distinguish the functions defined on the real number field M or 
the A^-dimensional Euclidean space M^, which is clear from the context. The L^(il inner product 
and the norm induced by the inner product are denoted by (•, •) and | • | respectively. The norm of 
the Sobolev space H"^{Q,) is denoted by || ■ \\m which is reduced to L'^{Q,) when m = 0. H^^{Q,) is 
the dual space of Hq{Q). We use D(0,T) to denote the C°°-functions with the compact on (0,T) 
and T>{0,Ty the space of distribution associated with P(0,T). 
Now we introduce some function spaces as follows. 



D 
H 

V - 



{i) E (C7~(f)))^| div(V'(3;) = on and 4>r{x) = on Fs}, 
completion of Dunder the {L'^{Q))^ — norm , 
completion of L'under the {H^{ft))'^ — norm , 



Dr, ={(pe C°°{n)\ ip{x) = on Fi}, 
H = closure of D^j in L^(il), 
W = closure of Dr, in H'^{n). 



(2.1) 



If there is no the anxiety 
• I and the norm of W is 



The norm of H is denoted by | • | and the norm of V is denoted by || • | 
confused in the notation, then we will also denote the norm of H is by 
denoted by || ■ ||. 

Suppose that {z,w) is classical solution of (jl.ip . 

Now, specially, we assume that viji G L'^{0;T; (L2(Fi))^) and vi,f2 G L'^{0;T; (^^(Fi))). 
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Multiplier the first equation of (jl.ip hy (j) ^ V, integrate by parts over and take the boundary 
condition into account to get 

^{z,ip) +a^(^)(z,V') + b{z,z,^) - {/3gw,<l)) = ifi,ip) + (i^i,V'n)ri, tpn = ii' ■ n)n. (2.2) 

Multiplier the second equation of (jl.ip hy ip & W , integrate by parts over F, and take the boundary 
condition into account to obtain 

^{w,ip) + ak^^'j{w,ip) + c{z,w,ip) = (/2,¥?) + (z^2,V')r2- (2-3) 

Now let X G Ci[0,T] be a function such that xiT) = 0. Multiplier Equations (H^]) and ([23]) by x 
respectively, and integrate by parts to yield 
rT 

-{zit),(j)x'it)) + a^(^^){z{t)Ax{t)) + Kz{t), z{t), <Px{t)) 



+{Pgw{t),<l)xmdt= r Ui{t).Mt))dt (2.4) 



L 



+ 11 {i^i,(pnXit))dsdt + {zo,(p)xiO). 
Jo JFi 

T 

[-{w{t), ^x'{t)) + ak(^)iw{t), (fixit)) + c{zit),wit), (f)xit)) 
]dt= [ {f2{t),^x{t))dt+ f [ {u2,px{t))dsdt + {wo,v)x{0), 



(2.5) 



'0 Jo Jr2 

where we write z{-,t) = z{t),w{-,t) = w{t) by abuse of notation without confusion from the context, 

and 

N 

{w{x, •)ip{x)dx, 



V / Zj{x,-)cl)j{x)dx, {w,ip)= / {i 
■ -^ Jn Jn 



0"y{w)iu,z) = {^{w)rotz{x, ■),Tot(l)), b{z,z,(p) = / (rotz(x, •) x z{x, ■))(j){x)dx, 



n 



ak{w){w,(p) = 22 / ^M^"^ a dx, c{z,w,(p) =22 / Zji^r) ^ ' ' (p{x)dx. 

j^i Jn ci^j -^^ Jn cixj 



We are now in a position to state some preliminary results. Throughout the paper, we always 
assume that there are positive constants 7o,7i,feo;^i such that 

70 < 7(0 < 71, A:o<A:(t) <A:i,Vte [0,oo). (2.6) 

The following Lemmas 12 . 1112 . 3 1 can be obtained by the Sobolev inequalities and the compactness 
theorem. We can also refer to theorem 1.1 of |14] on page 107 and lemmas 1.2, 1.3 of |14) on page 
109 (see also chapter 2 of |13]). The similar arguments can also be found in lemmas 1, 5 of 

Lemma 2.1. The bilinear forms a'^(.^.) and afc(.,.) are coercive over V and W respectivel That is, 
there exist constants c\,c\ > such that 



(•^)(z, z) > 7oCi||z|[f , V z and afc(^)(?i;, w) > /cqc'i V w G W. 
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Lemma 2.2. The trilinear form is a linear continuous functional with respect to each 

variable defined on [H^[Q))^ . That is, there exist a constant C2 > such that 

\b{u,v,w)\ < C2||n||||u||||tL'||, V u,v,w G {H^{Q,))^ . 

Moreover, the following properties hold true 
(i). b{u,v,v) = 0,y u,v e V ; 

(a). b{u, V, w) = —b{u, w,v),y u S V,v,w € {H^{ft))'^ ; 

(Hi). If Um u weakly on V and Vm ^ v strongly on H, then b{um,Vm,w) — )■ b{u,v,w),y w G 
V,weV. 

Lemma 2.3. The trilinear form c(-,-,-) is a linear continuous functional with respect to each 
variable defined on V x W x W . That is, there exist a constant C3 > such that 

\c{z, w,ip)\ < C3||z|| ||(/?||, V z G V,w,ip G W. 

Moreover, the following properties hold true 
(i). c{u,w,w) =Q,\I z ^V,w ; 
(a). c{z, w, if) = —c{z, ip, w), V z & V,w,ip € W ; 

(Hi). Ifzm — ^ z weakly on V and Wm — >■ w strongly on LP'{VL), then c{zm, Wm, ¥') b{z, w, f/?), V z G 
V,we L^{n),ip£ W. 

By considering (|2.3p . (|2.4p . we define the weak of (jl.ip such as; 

Definition 1. Let Y = Z x W = {L'^{0,T;V) n L°°(0, T; i/)) x {L'^{0,T;W) n L°^{0,T; H)). 
Suppose that 

/i G L^{0,T;V*),f2 G L^O,T;W*),vi G L^{0,T; {H-y^{T,))^),V2 G L^{0,T; {H~^I^{T2)), 

The pair {z, w} is said to be a weak solution of (|l.ip if it satisfies 
{z,w} £Y,z' e L\0,T;V*),w' G L\0,T;W*), 

{/,(/)) +a^(^)(z,0) + b{z,z,(j)) + {I3wg,(j)) = + (vi, 0n)ri , V; (/» G V, 

{w',(p) + afc(^)(?i;,v9) + c{z,w,(p) = {f2,^) + (^2, V3)r2, V; 99 G W, 
z{0) = 20,^(0) = Wo, 

where the inner product in is also denoted by (•, ■) without confusion from the context, 

and that in V* and V (also H~^{Q) and Hq{Q)) are denoted by (•, •). 

Next, we reformulate Equation (j2.7p in to the operator equation. To this purpose, it is noticed 
that for a fixed G F, the functional V) — )• a^(^)(z, (p) is linear continuous. So there exists an 
A^z G V* such that 

(A^z,^) = a^^.^){z,^),y (t>eV. (2.8) 
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Similarly, for fixed u,v € V, w{(z V) — ?• b{u,v,w) is a linear continuous functional on V. Hence 
there exist a B{u,v) G V* such that 



{B{u,v),w) =b{u,v,w),y w eV. (2.9) 
We denote B{u) = B{u,u). Define 

Lm = ifiA) + {viAn)r„y (peV. (2.10) 
Then Li a linear continuous functional on V and so there exists a constant C4 > such that 



11^1011 <c4<P\\v,y<p(^V. 

Hence there exists a Fi (^V* such that Li{(j)) = {Fi, (p) for all (j) . 

By the operators defined above, we can write the second equation of (12. 7p as 

dz 



(2.11) 



dt 



+ A^z + B{z) + (3gw = Fi. 



(2.12) 



Similarly, we have 



{AkW,if) = ak(^^'j{w,ip), {C{z,w),ip) = c{z,w,ip),AkW,C{z,w) e W*,y ip e W. (2.13) 



Define 



(2.14) 



Then L2 a linear continuous functional defined on W and so there exists a constant C5 > such 
that 

\\L2ip\\ < c^W^fWw ,y ^ & W. (2.15) 

Hence there exists a -F2 € W* such that ^2(95) = (-^2) V') for all <p € VF. By these operators defined 
above, we can write the first equation of (12. 7|) as 



— + AkW + C[z,w) = F2. 



(2.16) 



Combining (j2.12p and (|2.16p . we can write (j2.7p in the abstract evolution equation as follows: 



^ dz 

— + A^z + B{z) + Pgw = Fi, 

— + AkW + C{z,w) = F2, 

^ z{0) = zo,w{0) = Wq. 



(2.17) 



Lemma 2.4. If z G L'^{0,T;V), then B{z) G L^{{),T;V*); and if w G L'^{Q,T;W), then C{z,w) G 
L^{{),T;W*). 

Proof. Apply the Holder inequality and the Sobolev embedding theorem H^iVi) C L^(0) to obtain 

\{B{z),4,)\ = \h{z,z,ci,)\ < cq\\z\\4z\\l4Hl^ < C6lkll?||</'||i 
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and hence ||-B(z)||v* < ceH^Hf which shows that B(z) S L^{0,T;V*). Similarly, we have 

\{C{z,w),ip)\ = I - c{z,ip,w)\ < CY\\z\\L4\\ip\\\\w\\L4 < Cs\\z\\\\ip\\\w\\ < Cg{\\zf + || ^) || || 

and hence \\C{z,w)\\w* < cgdi-zlli + which shows that C{z,w) G L^{0,T; W*) for all tp €W. 

The proof is complete. □ 

We specified the constants Ci,i = 1,2, ••• ,9 are in this section in the remaining part of the 
paper. The following Lemma 12.51 comes from theorem 2.2 of [13j on page 220. 

Lemma 2.5. Let Xq,X,Xi be Hilbert spaces with the compact embedding relations 

Xq^ X ^ Xi. 

The for any bounded set K C R,v > 0, the embedding H'^{M.; Xq, Xi) C L^(IR;X) is compact, 
where 

H-^{R;Xo,Xi) = {ve L\R;Xo)\ snpp{v) C K,p^v G L\R;Xi)}, 

POO 

p'^v{t) = {2TTiTyv{T), v{t) = I e-2"*t;(t)dt, 

\\v\\h-^{R;Xo,Xi) 

3 Existence of the weak solution 

This section discusses the existence of the weak solution defined by Definition ([1]) to Equation (jl.ip . 
The main idea is to construct to a Galerkin approximation. 

Choose two orthogonal bases {uj}^^^ for V and for W respectively. Construct the 

Galerkin approximation solutions 

m m 
Zm{x,t) = '^qim{t)Uj{x), Wm{x,t) = himit) Hj (x) (3.1) 

i=i i=i 
such that for all j G {zm,Wm} satisfies 

{z'^{t),Uj) + a^(^^)(zm(t),nj) + b{zm{t),Zm{t),Uj) + (5{wm{t)g,Uj) = {fi{t),Uj) 

{w'^{t), flj) + a^,^^){Wm{t), fij) + c{Zm{t),Wm{t), fij) = (/2(t),/ij) + {v2,Hj)r2, 
Zm(0) = ZmO Zq in H, Wm{0) = WmO Wq in H J = 1,2 ■ ■ ■ ,m. 

where z^o is, for example, the orthogonal projection in H of zq on the space spanned by ui, ti2, • • • , Um 
and WmO is the orthogonal projection in H of on the space spanned by /^i, /i2, • ' ' , ^^m■ 

Hence once again, we write Zm{t) = Zm{-,t),Wm{t) = Wm{-,t), fi{t) = fi{-,t),Vi{t) = Vi{-,t),i = 
1, 2 by abuse of notation without the confusion from the context. It is seen that for any m G N"*", sys- 
tem ()3.2p is a system of nonlinear differential equations with the unknown variables {qjm{t)^ hjm{t)} 
and the initial values qjmiO) = izo,Uj), hjm{0) = (wq, fJ-j), j = 1, 2, • • • , m. By the assumption, this 
initial value problem admits a solution in some interval [0,tm]- We need to show that t^ = T. 
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Lemma 3.1. Let {zm,Wm} be the sequence satisfying 13. Then there exists a subsequence of 
{zm,Wrn}, still denoted by itself without confusion, such that 

Zm ^ z weakly in L^(0, T; V) and z^ ^ z weakly star in L°°(0, oo; H), (3-3) 

where z e L^iO,T;V)nL°°{0,T;H), and 

Wm — > w weakly in L^(0, T; W) and Wm w weakly star in L°°(0, oo; H), (3.4) 

where w € ^^(O, T; VF) n L°°(0, T; i?). 

Proof. By Lemmas 12.2112.31 we have 

{z'j^{t),Zm{t)) + a^(^yj^){Zm{t), Zm{t)) + f3{Wm{t)g, Z.m{t)) = {fi{t),Zm{t)) 

(t))r,, (3.5) 
By assumption ()2.6p . for any given e > 0, we can get from (j3.5p that 



j^\Zm{t)\'^ + 2jo\\Zmit)f < -2/3{Wm{t)g,Zmit))+2{fi{t),Z„,{t)) 



+2{viit),z^^{t))r, <7olk. 



+ :^ll/i(i)lly*+/3||5lloo 0lkm(t)f + 
+£^||^m„(i)II^i/2(r^) + ~2lbi(0ll//-i/2(ri)- 



(3.6) 



Here and hereafter, ||2;n"i(t)||^^i/2(ri)]iv we denote by |knm(i)||^i/2(r^] and ||fi(t)||^^_i/2(r^)]jv by 
lbi(*)ll^-i/2(rj) simply. 

By the trace theorem from H^{il,) to H^/'^{T), there exists a constant cio > such that 



1 1 ^ITln 

By substituting above inequahty into ()3.6p . we obtain such as; 



Jt\^m{t)\' + [70 - mWooe' + cioe2)]||z„(t)f 

< ^\\fi{t)fv* +m\oo^\Mt)f + ^|bi(t)||^-v2(ro- 

To ^ ^ 

Setting £2 = 7o/(2(/3||g||oo + cio) in ([32]) gives 

^Mt)? + f lk™(t)f < + '^^"'"°°^'^°V l|g||oo||||t.^(t)||^ 

2 7oCi 7oCi 

, 2(/3||g||oo +cio) 2 

By assumption (12. 6p again, for any given e > 0, we can get from ()3.5p that 
^|ti;^(t)|2 + 2/co||w;.„(t)f < A;o||t«„,(t)f + ■^||/2(t)||vF* 

-2|i ., /.M|2 , 1 11 . /.M|2 



(3.7) 



(3.8) 



(3.9) 



I|w^m(t)||j;^l/2(r2) + -2\\v2{t)\\H-^/2{T2)- 



By the trace theorem from H^{il.) to H^^^{T) again, there exists a constant cn > such that 

(0lli/l/2(r2) ^ Cii\\'Wm{t)\\. 

By substituting above inequality into (13. 9p . we obtain such as; 

+ (fco - c^,e')\\wm{t)f < Ijmml.. + ^Ib2(i)||^-v2(r,)- (3-10) 
Setting £2 = A;o/(2cii) in (IXTIl gives 

+ y < pl!/2(t)fvi^* + ^Ib2(t)||^-v2(r,)- (3-11) 



..- .loo ,, 

70 Jo 



70 Jo 
Therefore 

'2 



eS.SUp|z„(t)|2 <- r||/l(t)f^.cit+ ^^^"^"°°+"^°V ||g|U|| / \\Wrnit)fdt 



t 70 Jo 70 Jo 

+*i^rii«(t)iii,...„.,^*+NP 



(3.12) 



Integrate (|3.1ip over [0, T] with respect to t to give 

+ |^UO)P<f r\\f2{t)fw,dt + ^ rik(i)||^-V2(r,)rft+K0)|2. 

i/ i/ 

Since the right-hand side of ()3.12p is bounded, we have 

{wm} is a bounded sequence in L^(0, T; 1^). (3.13) 
Replace T by t e [0, T] in (f3l2]) to obtain 

ess sup \wm{t)\^ <^ [ \\f2mw*dt+'^ [ \\v2ml.^;,.^.dt + \w{0)\\ (3.14) 
t «^o Jo £ Jo 

Hence 

is a bounded sequence in L°°(0, T; L^(r2)). (3.15) 
On the other hand, integrate (j3.8p over [0, T] with respect to t to give 

km(T)p+5 r \\z^{t)fdt <- r \\h{t)\\i,dt 

^ Jo 70 Jo 

+ " ^^"u;^(Ofcii (3.16) 



{zm} is a bounded sequence in L (0,T;y). (3-17) 
Replace T by t e [0, T] in (ITOD to get 



(3.18) 



70 

Therefore, 

{zm} is a bounded sequence in L°°{0,T; H). (3.19) 
iK3\i and ([331) then follow from (I3T3I1 , (ISTTll (l3T5]l and (lOTIl . □ 



Lemma 3.2. Let {zm,Wm} be the sequence determined by Lemma \3.1[ Then there exists a sequence 
of{z rajWrn}; stUl denoted by itself without confusion, such that 

Zm ^ z strongly in L^{0,T; H), Wm w strongly in L'^(0,T; H). (3.20) 

Proof. By virtue of Lemmas 12.1112.31 we can write p.2p as follows: 

/ dZm{t) 



dt 



{Fi - I3gwm{t) - A^Zm{t) - B{zm{t),Uj),\f j = 1,2,--- ,m, 



^^"dt^ = (-^2 - C{Zm{t),Wm{t)) - AkWm{t) , flj) ,\/ j = 1,2, • • • , m. 



(3.21) 



Denote by {zm-,Wm} the {zm,Wm} with zero values outside of [0,T] and {zm,Wm} the Fourier 
transformations of {zrmWm}- We claim that there exists a > such that 

oo 

\t\'^'' \\Zm{T)f dT < oo. (3.22) 

-oo 

To this end, we write the first equation of (I3.2ip as 

^{Zm,Uj) = {flm,Uj) + {zom,Uj)5o - {Zm{T) ,Uj)5T , (3.23) 

where 5q,5t are Dirac functions, and 

fim{t) = fimit) for t G [0,r] and /i^(t) for t > T, 

flm{t) = Fi- (3gWm{t) - A-yZm{t) - B{Zm{t)). 

Take Fourier transform for Equation ()3.23p to get 

2TriT{zm,Uj) = (/im, Uj) + {zom, Uj) - (^^(r), Uj)e"^"^'^, (3.24) 

where fi fn is the Fourier transform of f\m- 

Let qjmit) be the function of qim in Zm{t) = X^jli 1jm{t)uj that is zero outside of [0, T] and let 
Qjmi'T) be its Fourier transform. Multiplier Equation (|3.24|) by qjm and sum for j from 1 to m to 
obtain 

27ri\zm{r)\^ = {fim{r),Zm{T)) + (zom(0),5„^(r)) - (z^ (T) , i„ (t) )e-2-^\ (3.25) 
We thus conclude that 

flm{t)\\v''dt < I [||/l(t)||v* +C2||u'^(t)|| +7l|kmWII +Cl||Zm(t)f] dt, (3.26) 

Jo 

where we used the fact that ||i?2:m(i)|| < ci ||2;m(t) |P- By (j3.13p . (|3.15p and (j3.17p . it follows from 
(|3:26D that 

sup||/i^(r)||y* <oo. (3.27) 
Apply ()3.27p and the facts sup„g^+ [||zm(0) || + ||zm(r)|l] < oo to (|3.25p to yield 

\T\%n{T)\^ < C3pm(T)|| + C4|5m(r)| = C5 ||zm(r) || , C5 = c'3 + C4. (3.28) 
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For V , fixed < 1/4 , we can observe that 
From this inequaUty, we obtain 



oo 



T\^''\zm{T)\^dT<c'{u) I -^±^\z^{r)\^ dv (3.29) 



1 + r, 

-oo J — oo ^ 1 ' I 

Thus by ()3.28p . we can obtain the constants Cq > and c'l > such that 



oo roo 1 + iT- 

ool+kl 

J —oo ^ ' Pi J —oo 

By (I3.20|) and the Parseval equahty the last integral is bounded as m — >■ oo ; thus p.22p will be 
proved if we show that; 

\^m{r)\\ Const 



1 + |t|1-2'^ 

By the Schwarz inequality and the Parseval equality we can obtain 

which is finite since < 1/4 , and bounded as m — )• oo ; by (13.19p . The proof of p.22p is achieved. 
By ([3:22]) . (ISTTll and (l3T9]l . we conclude that 

{z„,} is bounded in H^iR; V) n i7''(M; F). (3.30) 

By Lemma 12.51 there exists a subsequence of {zm,Wm} that is still denoted by itself without 
confusion such that 

Zm^ z strongly in L'^{0,T;H), Wm ^ w strongly in L'^{0,T; L'^{Q)). 

This is ([3:20]) . □ 

Theorem 3.1. Suppose that the functions 7, A; satisfy that when Wm w in L'^{0,T; L^{Q)), then 
li'Wm) ~^ l{w).,k{wm) — > k{m) in L'^{0,T; L'^(Q,)), respectively. Then there exists a weak solution 
to / TO) . 

Proof. Let ^ and 9 be continuous differentiable vector functions defined on [0,T] with ^(T) = 
9{T) = 0. Multiply the first equation of ()3.2p by ^ and integrate over [0, T] with respect to t to 
give 



{Zm{t), '^'{t)uj)dt + / [a^(^^)(z„(t), '^{t)Uj) + b(z.m(t),Z.m{t), '^{t)uj) 

^ (3-31) 

+(3{w^{t)g,^{t)uj)]dt = {zQm,Uj)^{Q)+ / {fiit),uj)dt+ / {vi,^{t)ujjridt. 
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Multiply the second equation of (j3.2p by 9 and integrate over [0, T] with respect to t to give 
rT r-T 

{Wmit),9'it)flj)dt+ / [ak(wm)i'^m{t),0it)fij) + c{z,nit),Wmit),9{t)fIj) 

^ rr. (3.32) 

]dt = {wo^,fXj)9iO)+ {f2{t),9it)fij)dt+ {v2,9it)fij)T,dt 
Jo Jo 

Passing to the limit as m — > oo in (|3.3ip and (j3.32p by applying (|3.3p . (j3.20p . the properties (iii) 
in Lemmas 12.21 and 12.31 for b and c, and the continuous assumption on 7 and k, we obtain 



{z{t), ^'{t)uj)dt + I [a^^^^){z{t),^{t)uj)dt + h{z{t),z{t), ^{t)uj) 
Jo 

+P{w{t)g,-^{t)uj)]dt = {zQ,Uj)-^{^)+ [ {h{t),uj)dt+ [ {vi,^it)ujjr,dt. 

Jo Jo 

{w{t),9'{t)iij)dt+ [ [ak(^^)iwit),9it)fij) + ciz{t),wit),9it)fij)]dt 
Jo 

= {wo,fij)9{0)+ [ {f2{t),9{t)fij)dt+ [ {V2,9{t)fij)r,dt, 
Jo Jo 

where in obtaining (|3.33p and (j3.34p . we used the following facts: 



(3.33) 



(3.34) 



The convergence of the nonlinear terms in b{-, •, •) and c(-, •, •) can be obtained in the same 
way as that in Chapter 3 of [13J. 



0"riwm)i^rn{t),'^{t)fij)dt = / {j{Wjn)rotZm{t),'^{t)rotUj)dt 



Jo 



T 

{rotZm{t),^{wm)i'ot^{t)rotUj)dt 





T r-T 

{rotz{t),'y{w)^{t)rotUj)dt = / {'y{w)rotz{t),'^{t)rotUj)dt 
^0 







T 

a'r{w){z{t),^{t)fij)dt, 



where we used the facts that rotzm rotz weakly in L (0,T;y) and j{wm)^it)rotu 



'J 



'y{w)^(t)Vuj strongly in L^{0,T; H) from the assumptions of Theorem 13.11 
• Similarly 

/ ak{w„,){wm{t),9{t)^j)dt^ I a^(^^){w{t),9{t)fij)dt 
Jo Jo 

by the facts again Vwm Vw weakly in L^{0,T;W) and k{wm)9it)V fij k{w)9(t)Viij 
strongly in L^(0,T; H) from the assumptions of Theorem 13.11 

By the density arguments, we have that (|3.33p and (|3.34p hold true for any cj) € V instead of 
Uj and if £ W instead on fij, respectively. That is, 

iz{t),^'{t)(l))dt+ [ [a^(^^)izit),^it)cl))+b{z{t),zit),^it)cl)) 

(3.35) 



JO 



+(3{w{t)g,^{t)cP)]dt = (zo, 0)^(0) + r {h{t),4>)dt + r{vi,^{t)4>n)r,dt,y4>e V- 

Jo Jo 
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{w{t),9'{t)^)dt + 1 [afe(^)(u;(t), e{t)^) + c{z{t),w{t),e{t)^)]dt 

° ^ ^ (3.36) 

= {wo,^)9{0)+ {f2it),9{t)ip)dt+ {v2,e{t)ip)r,dt,y^GW. 
Jo Jo 

Now take G {D{0,T))^ in (I335D and G l?(0,r) in (l336]l . Then {z,^} satisfies 

f (2',<^) + ^(^'^'<?^) + i(^wg,(t>) = if 1,(1)) + (i;i,(/)n)ri,V (/> G 

1 (w^ip) + afc(^)(iz;,((3) + c(z,u;,v3) = {f2,^p) + (^2, 99)r2, V G W. 

This is the equations in (j2.7p . Finally, we determine the initial value of Actually, multiply 

the first equation of (j3.37p and integrate over [0, T] with respect to t to get 

zit),^'it)(f>)dt+ 1 [a^(^)(z(t),^(t)0)(it + 6(z(t),z(t),^'(t)(/)) 

^ ^ ^ (3.38) 

+P{w{t)g,^{t)<t>)]dt = {z{d)A)^{Q) + / {h{t)A)dt+ / (7;i,M/(t)</)„)ri(it,V(/)Gy. 

Jo Jo 

Subtract ([3381) from (l335]l to get (z(0) - zo,'/')^(0) = 0. Take ^' so that ^'(0) = 1 we get 
{z{Q) — zq, (J)) = for all (p & V. So z(0) = zq. The similar arguments lead to w{0) = wq. The proof 
is complete. □ 



4 Uniqueness of the weak solution 

Theorem 4.1. Let 7, A; : L^(r2) — ?> L°°(il) are Lipschitz continuous functions, that is, there are 
constants /ii,/i2 > such that 

Il7('"^*) — 7(if**)||L°° < ^ill'tL'*) — w^,^\\,t G [0,T] a.e., V ti;*,?/;** G L'^{^), 

(4.1) 

||A;(t(;^,) — A;(ti;*,,)||Loo < /2|[w;*) — tf^**||,t G [0,r] a.e., V if*, w** G L^(ri). 
Suppose that the weak solution {z,w} of claimed by Theorem \3.1\ satisfies 

Ad(^^^^^ + -^\\w{t)\\l,) <7o,VtG[0,T] a.e., (4.2) 

where d > is the constant in Sobolev inequality that \\f\\L*(^n) — for all f G H^{i}), 

which depends on N and fl, and the constants ci and c\ are that in Lemma \2.1\ Then the weak 
solution is unique. 

Proof. Suppose that we have two weak solutions {z^,w^,}, {z^,,z^^} to (jl.ip . Set z = — z^^,,w = 
Tu* — tu**. Then by (12. 7|] 

' d 

— iz,<j)) +{'y{w^)Vz, + b{z^,z^,(t)) - b{z^^,z^^:,<j)) + {I3wg,4>) 
dt 

= {{l{w^) - 7('u;**))Vz*^,, V0),V G y, 

d ^^-^^ 

— {w,(f) +{k{w^)Vw,V(p) + c{z^,w^,,(f) - c{z^^.,w^.^,(p) 

= —{{k{w^) — A;(w**))V?i;**, V99), V tp e W. 
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Set cp = z,(p = w in ()4.3p to get 

r d 



dt 



d 



(z, z) +{'-^{w^)V z, Vz) + b{z^,z^,z) - b{z^,^,z^,^,z) + {fiwg, z) 
= {{{l{w*) - 7(w**))V2:**, V(?!)), 



(4.4) 



= — — A;(w**))V?i;**, Vw). 
This together with b{z**,z,z) = c(z,,*, tf) = gives 

\z{t)f + 27o|lVz(t)f < 2b{z, z,z^) + 2{/3gw, z) - 2{{-i{w^) - 7(u;„))Vz„, Vz), 
\w{t)f + 2ko\\Vw{t)f < 2c{z,w,w^) + 2{{k{w:,) - k{w^^))Vw^^,Vw). 
By the Lipschitz continuity (j4.ip . for any given e > 0, it follows from the "z part" of ()4.5p that 

'2 + 27o||Vz(t)f < 2||z||i4||Vz||||z*||i4 + 2(3\\g\\ 



— I 
eft' 



(4.5) 



dt' 

-7(^t'**)l|L°°||Vz„||||V2;|| 

< 2||z||2,4||Vz||||z*||j;^4 + 2/3||5(||oo||w||||-z|| + 2/i||?i;* - t(;^,*||i2 ||Vz** || ||Vz| 



(4.6) 



< 2||z||i4||Vz||||z,||i4 +2/3||5||oo||u;||||z|| +2/i 
Putting = ^2//^ in (ITH]) . we obtain 

^||z(t)f + 27o||Vz(t)f < 2||z|[i4||Vz||||z,||i4 + 
Similar arguments to the part" of (j4.5p . we have 



70 



^Ili^llllVz^f + i^||Vz 



70 



Ui^ \\z\\ + 



270 ci 



|w'lli2||Vz„||. (4.7) 



dt 



|u;(t)r + 2A;o|lVw(t)r < 2||z||i4 1| Vu;|l Hu-* 11^4 + 



2koCi 



\w\\l2\\Vw 



(4.8) 



Sum (I12D and (gSD to get 



d 



j^iimf + \\w{t)f]+jo\\S/z{t)f + ko\\Vw{t)f 



< 2\\z\\l4\/z\\\\z4h + 



w z 



+ o ^ II^^IIl^II^^**!! + 2|k||L4||Vu;||||t(;*||i4 
^70 ci 



(4.9) 



Apply the Holder inequality and the Sobolev inequality to the right-hand side of (|4.9p to obtain 
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for any given ei > 0, that 



< 2d|iVzf ||z*||i4 + 2/3||ff||oo||u;||||2|| + -^||t/;||i2||Vz„|| +2d||V^||||Vw;||||u;*||i4 

270 Cl 



+^^II^IIl2||Vu;„|| (4.10) 

< 2d|iVzf ||z,||i4 + 2/3||g||oo||u;||||-z|| + ^T^lMliWVz^^l 

270 Cl 

+2d(^^\\Vwf + ^J\Vzf\\w4h^ + 



^2 II ||2 IIV7 Il2 



2kQCi 

where d > is the constant from Sobolev inequality that ||/||L4(f^) < j:/i(f^) for all / G H^{Q,). 
Set ef = koCi/2d in ()i30]) to get 

+ + 7oci||.(i)f + ^IKt)f 

< 2d\\zf\\z4L4 + p\\g\\oo{\w\'^ + \z\^) + TT^kPlk**!! + T^\\zf\\\w*\\l4 

270C1 fcoci 



I '2 I |2|| |t2 

+ , \w\ \\w^^\\ 
2feoCi 

2d\\z*\\L4 + -rr^\\w*\\'i4 ) \\z\\'^ + ( /3||5(||oo + 2 11^**11^ + , ) 



(4.11) 



+/3||5llookP 

^ 70Cl II ||2 , ( n\\ II , ^1 II ||2 , h It ||2\ I |2 , 

2 V 7oCi 2feoCi y 

In the first term of the last row of (|4.1ip . we used the assumption (|4.2p . From (j4.1ip . we conclude 
that 

|[Kt)p + KOp] + ^||z(t)f + ^IHt)f 

V 270C1 2fcoCi ) 



from which we obtain 



j^\z(t)f + < M(t)|u;(t)p + N\z(t)f (4.12) 



or 

|[|z(t)P + |^(t)p] < (M(t) + N)\\w{t)f + |z(t)p), (4.13) 

where 

M{t) = m\oo + i^\Mt)f + -^\\w,4t)f,N = P\\g\\^. 
270C1 2feoCi 

Since M{-) + is integrable in [0,r] with respect to t, we obtain, from ()4.13p . that 

|{e-(^'^+^)[|z(t)|2 + |t,(t)|2]} <0,Vt G [0,r] a.e.. (4.14) 
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This together with z(0) 



w(0) = gives 



+ \w{t)\^ < 0,Vt G [0,r] a.e.. 



Therefore, z^. 



= w^^:. The proof is complete. 



□ 



Remark 4.1. If we denote (see, e.g., [IT]) 



Re(t) 



112,4 (the Reynold number); Ra(t) 




(the eigh number) 



7oCi 



Then, condition ()4.2p is reduced to condition: 



Re{t) + Ra{t) < 1 



Remark 4.2. Theorem 14. II is a generalization of the results of [T7]. In [T7j, the boundary condition 
for the velocity of fluid is given by the standard boundary condition, that is, the homogeneous 
Dirichlet boundary condition, and the boundary condition for the temperature of fluid is given by 
the non-homogeneous Dirichlet boundary condition. 
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